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Recent progress in the determination of \V US \ employing strange hadronic r-decay data are reported. This 
includes using the recent OPAL update of the strange spectral function, as well as augmenting the dimension-two 
perturbative contribution with the recently calculated order al term on the theory side. These updates result 
in \V US \ = 0.2220 ± 0.0033, with the uncertainty presently being dominated by experiment, and already being 
competitive with the standard extraction from K B z decays and other new proposals to determine \V US \. In view 
of the ongoing work to analyse r-decay data at the B-factories BaBar and Belle, as well as future results from 
BESIII, the error on \V US \ from r decays is expected to be much reduced in the near future. 
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1. INTRODUCTION 

In the past decade hadronic r decays have been 
an extremely fruitful laboratory for the study of 
low-energy QCD. Detailed investigations of the r 
hadronic width 
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as well as invariant mass distributions, have 
served to determine the QCD coupling a s to a 
precision competitive with the current world av- 
erage [1-5]. The experimental separation of the 
Cabibbo-allowed decays and Cabibbo-suppressed 
modes into strange particles [6-10] also opened a 
means to determine the mass of the strange quark 
[11-17], one of the fundamental QCD parameters 
within the Standard Model. 
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These determinations suffer from large QCD 
corrections to the contributions of scalar and 
pseudoscalar correlation functions [1,11,18,19] 
which are additionally amplified by the partic- 
ular weight functions which appear in the r 
sum rule. A natural remedy to circumvent 
this problem is to replace the QCD expressions 
of scalar and pseudoscalar correlators by corre- 
sponding phenomenological hadronic parametri- 
sations [6,13,15,20,17], which turn out to be more 
precise than their QCD counterparts, since the 
dominant contribution stems from the well known 
kaon pole. 

Additional suppressed contributions to the 
pseudoscalar correlators come from the pion pole 
as well as higher exited pseudoscalar states whose 
parameters have recently been estimated [21]. 
The remaining strangeness-changing scalar spec- 
tral function has been extracted from a study of 
S-wave Kir scattering [22,23] in the framework of 
resonance chiral perturbation theory [25]. The 
resulting scalar spectral function was also em- 
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ployed to directly determine m s from a purely 
scalar QCD sum rule [24]. 

Nevertheless, as was already realised in the first 
works on strange mass determinations from the 
Cabibbo-suppressed r decays, m s turns out to de- 
pend sensitively on the element | V us | of the quark- 
mixing (CKM) matrix. With the theoretical im- 
provements in the r sum rule mentioned above, 
in fact \V US \ represents one of the dominant un- 
certainties for m s . Thus it appears natural to 
actually determine | V us | with an input for m s as 
obtained from other sources [17,26,27]. 

Succeeding the high-precision status on r- 
decay observables already attained by by ALEPH 
and OPAL at LEP and CLEO at CESR, now the 
B-factories BaBar and Belle are starting to pro- 
duce their first results on hadronic r decays, and 
in particular on Cabibbo-suppressed modes [28- 
30] . These two facts make the strange hadronic r 
decay data an ideal place for determining SU(3) 
breaking parameters such as | | and/or m s . 
The obvious advantage of this procedure is that 
the experimental uncertainty will eventually be 
reduced at the B-factories and at future facilities 
like the r-charm factory BEPCII 

2. THEORETICAL FRAMEWORK 

Employing the analytic properties of two-point 
correlation functions for vector [J = V) and 
axial-vector ( J = A) two quark-currents, 
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one can express R T as a contour integral running 
counter-clockwise around the circle |s| = M 2 in 
the complex s-plane: 
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Here, we have used integration by parts to rewrite 
R T in terms of the logarithmic derivatives 
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Moreover, experimentally R T can be decomposed 
into the following three contributions 

Rr = Rr,V + Rr,A + Rr,S , (5) 

according to the quark content 

IL 7 ( S ) = \V ud \ 2 {U J v ^ d ( S )+Ui ud ( S )} 

+ \Vus\ 2 {Tl J v ^( S )+Kus(s)} , (6) 

where R T y an d R t ,a correspond to the first two 
terms in the first line and R T ,s to the second line, 
respectively. 

Additional information can be inferred from 
the measured invariant-mass distribution of the 
final hadrons, which defines the moments 

At large enough Euclidean Q 2 = — s, both 
n i+T (Q 2 ) and T1 L (Q 2 ) can be organised in a di- 
mensional operator series using well established 
QCD operator product expansion (OPE) tech- 
niques. One then obtains 

R? = N c S EW {(\V ud \ 2 + \V us \ 2 ) [l + 6 kl ^] 

|K S | 2 « D) ]}.(8) 
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The electroweak radiative correction Sew = 
1.0201 ± 0.0003 [31] has been pulled out ex- 
plicitly and S kl ^ denote the purely perturba- 
tive dimension-zero contributions. The symbols 



kl(D) 



stand for higher dimensional corrections 



in the OPE from dimension D > 2 operators, 
which contain implicit suppression factors 

[1,11,13,18]. The most important being the oper- 
ators m 2 with D = 2 and m s (qq) with D = 4. 

In addition, the flavour-SU(3) breaking quan- 
tity 
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enhances the sensitivity to the strange quark 

kl(2) 

mass. The dimension-two correction 5- • is 
known to order for both correlators, J = L 
and J = L + T [11,13,32]. 

In Ref. [11], an extensive analysis of this D = 2 
correction was performed and it was shown that 
the perturbative J — L correlator behaves very 
badly. The J = L+T correlator was also analysed 
there to order a 2 and showed a relatively good 
convergence. In the following, we have included 
the recently calculated 0{a? s ) correction for J = 
L + T [32] . One can see that the J — L + T series 
starts to show its asymptotic character at this 
order, though it is still much better behaved than 
the J — L component. Due to the asymptotic 
behaviour, it does not make much sense to sum 
all known orders of the series. This question will 
be investigated in more detail by us in the future. 

3. DETERMINATION OF \V ua \ WITH 
FIXED m s 

One can now use the relation 
R m 

\V us \ 2 = Wo ^ , (10) 

IT/ 12 ott T,th 
\Vud\ 

and analogous relations for other moments, 
to determine the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element \V US \. Notice that on 
the right-hand side of (10) the only theoreti- 
cal input is 5R®° th , which is around 0.24 and 
should be compared to the experimental quan- 
tity R T ° v+A /\V ud \ 2 which is around 3.7. There- 
fore, with a not so precise theoretical prediction 
for 5R®° th one can get a quite accurate value for 
\V US \, depending on the uncertainty provided by 
experiment. 

The very bad QCD behaviour of the J = L 
component in 6R^\ h induces a large theoreti- 
cal uncertainty, which can be reduced consid- 
erably using phenomenology for the scalar and 
pseudo-scalar correlators [17,26,27]. In particu- 
lar, the pseudo-scalar spectral functions are dom- 
inated by far by the well-known kaon pole, to 
which we add suppressed contributions from the 
pion pole, as well as higher excited pseudo-scalar 
states whose parameters have been estimated in 



Ref. [21]. For the strange scalar spectral function, 
we employ the result [22] , obtained from a study 
of S-wave Kir scattering within resonance chiral 
perturbation theory [25] , which has been recently 
updated in Ref. [33]. 

The smallest theoretical uncertainty arises for 
the kl = 00 moment, for which we get 

SR°° th = 0.1544 (37) + 9.3 (3.4) m 2 s 

+ 0.0034(28) = 0.240(32), (11) 

where m s denotes the strange quark mass in units 
of GeV, and in the MS scheme at a renormalisa- 
tion scale of n = 2 GeV. The first term contains 
the phenomcnological scalar and pseudo-scalar 
contributions, the second term contains the rest 
of the perturbative D = 2 contribution, while the 
last term stands for the rest of the contributions. 
Notice that the phenomenological contribution is 
more than 64% of the total, while the rest comes 
almost from the perturbative D = 2 contribu- 
tion. Here, we update 5B™ th of Refs. [17,26,27] 
in various respects. Firstly, we use the recently 
updated scalar spectral function [33]; secondly, 
we include the corrections to the J = L + T 
correlator as calculated in Ref. [32]; and finally, 
we use an average of contour improved [34] and 
fixed order perturbation results for the asymptot- 
ically summed series, in order to have a more con- 
servative account of uncertainties resulting from 
unknown higher orders. A detailed discussion of 
this contribution will be presented elsewhere [35] . 

For the m s input value, we use the recent av- 
erage m s (2GeV) = (94 ± 6) MeV [33], which in- 
cludes the most recent determinations of m s from 
QCD sum rules and lattice QCD. The strange 
quark mass uncertainty corresponds to the most 
precise determination from the lattice. 

Recently, Maltman and Wolfe have criticised 
the theory error we previously employed for the 
D = 2 OPE coefficient [38]. Awaiting a more de- 
tailed study [35], in our updated estimate (11), 
we have decided to include a more conservative 
estimate of unknown higher-order corrections by 
using an average of contour improved and fixed- 
order perturbation theory. Still, we do not think 
that artificially doubling the perturbative uncer- 
tainty, as was done in [38], represents an error 
estimate which is better founded. Notice further- 
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more, that SR®° th is dominated by the scalar and 
pseudoscalar contributions which are rather well 
known from phenomenology, and that the larger 
perturbative uncertainty is compensated by the 
smaller m s error, so that our final theoretical un- 
certainty is practically the same as in previous 
works [17,26,27]. 

In order to finally determine \V US \, we employ 
the following updates of the remaining input pa- 
rameters: \V ud \ = 0.97377 ± 0.00027 [36], the 
non-strange branching fraction R^v+a = 3.471 ± 
0.011, ,[10] as well as the strange branching frac- 
tion = 0.1686±0.0047 [10] (see also Refs. [6] 
and [8]), which includes the theoretical prediction 
for the decay B[t -» Kv T {>y)] = 0.715 ± 0.003 
which is based on the better known K — > ^^(7) 
decay rate. For \V US \, we then obtain 

\V US \ = 0.2220 ±0.0031 cxp ± 0.0011th- (12) 

The experimental uncertainty includes a small 
component from the error in | V u d | , but it is dom- 
inated by the uncertainty in R^P S , while the the- 
oretical error is dominated by the uncertainty in 
the perturbative expansion of the D = 2 contri- 
bution. 

4. SIMULTANEOUS FIT OF \V UB \ AND 
m s 

In principle, it is also possible to perform a si- 
multaneous fit to I V us I and m s from a certain set 
of {k, I) moments. As soon as more precise data 
are available, this will be the ultimate approach 
to determine \V US \ and m s from hadronic r de- 
cays. With the current uncertainties in the data 
and a persistent question about a monotonous k- 
dependence of m s [17,27], a bias could be present 
in the method. Furthermore, the correlations be- 
tween different moments are rather strong and 
also have to be properly included on the theory 
side. 

Here, we shall restrict ourselves to a simplified 
approach where all correlations are neglected. For 
the simultaneous fit of |V^ S | and m s , we employ 
the five Rf moments (0, 0) to (4, 0) which have 
also been used in our previous analyses [17,27]. 
Performing this exercise, for the central values 



we find: 

\V US \ =0.2196, m s (2GeV) = 76MeV. (13) 

The expected uncertainties on these results 
should be smaller than the individual error given 
in eq. (12) and the one for m s presented in 
Ref. [27], but only slightly since the correlations 
between different moments are rather strong. 

The general trend of the fit result can be under- 
stood easily. m s from the simultaneous fit turned 
out lower than the global average m s (2 GeV) = 
94 ± 6MeV considered above. Thus, also the 
corresponding 5R T , t h is lower, resulting in a re- 
duction of |V^s|. Furthermore, the moment- 
dependence of m s is reduced as compared to our 
analysis [17] on the basis of the ALEPH data 
alone. Nevertheless, we shall leave a detailed er- 
ror analysis for a future publication. 

5. CONCLUSIONS 

High precision Cabibbo-suppressed hadronic r 
data from ALEPH and OPAL at LEP and CLEO 
at CESR already provide a competitive result for 
\V US \- As presented above and in Refs. [17,26,27], 
the final uncertainty in the t determination of 
\V US \ becomes an experimental issue and will 
eventually be much reduced with the new B- 
factory data [28-30], and further reduced at fu- 
ture r facilities. A combined fit to determine both 
\V US I and m s will then be possible. Hadronic r de- 
cays have the potential to provide the most accu- 
rate measurement of | V us | and a very competitive 
m s determination. 
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